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Low dimensional strongly correlated electron systems are of considerable recent interest 
as possible models of high-T c superconductivity. One dimensional versions of such models 
have proven to be especially interesting as they are often completely integrable and may 
be solved exactly. A particular case is the one dimensional Hubbard model, which was 
solved by Lieb and Wu in 1967 ^@|, and which exhibits a hidden 5*0(4) symmetry HIJ. 
Recently, an extended Hubbard model was introduced which is solvable by Bethe ansatz in 
one dimension |5]||. in this letter, we discuss the long range version of the one dimensional 
extended Hubbard model on a uniform lattice. 

There has been renewed interest in the integrable electronic models of Calogero- 
Sutherland type of strong correlations since the independent work of Haldane and Shastry 
IIU on the one dimensional spin chain of 1/r 2 exchange interaction flO|-19 . The system 



we consider here is the Jastrow-wavefunction-integrable extended Hubbard model defined 
on a one dimensional closed uniform chain of size L. This model is a generalization of the 
Haldane- Shastry spin chain and the long range supersymmetric t-J model and is similarly 
completely integrable. In the following, we will investigate the ground state and excitations 
of the system in detail. 

In the extended Hubbard model, each lattice site may be either unoccupied (hole), oc- 
cupied by a single electron (up-spin or down-spin) or doubly occupied. The latter case of 
localized electron pairs are denoted localons. Thus there are four possible states at each site 
— two bosonic and two fermionic. The Hamiltonian then takes the form 

H = ~\ E J i^ij, (1) 

where the permutation operator exchanges the states between sites i and j. Uij is given 
explicitly in terms of electron operators by 

% = c jT c it( 1 _ n a ~ n 3i) + c it c jT( 1 - n n - n n) 
+c n c n( 1 ~ n ^ ~ ra ;t) + c h c ji( l ~ n ^ - n it) 

+ -(m - l)(nj - 1) + cl^Cj^ + c ll c i1 c] 1 c] l 



--(n iT -nn)(n j} -n si ) 
~ c h c i1 c h c jl ~ c h c il c ll c tf 

+ (n ir - l/2)(n a - 1/2) + (n jT - l/2)(n, u - 1/2). (2) 

The operators c icr and cj CT are electron annihilation and creation operators at site i with 
spin a =|, j. The electron number operators are given by = cLc^, = cjiqi and 

Because of the fermionic nature of the electrons, picks up an appropriate sign when 
exchanging fermions, and hence acts as a graded permutation operator. This leads to an 
S77(2|2) supersymmetry of the system ||. In fact there is a natural generalization of this 
model to a chain where there are m bosonic and n fermionic species at each site, leading 
to an SU(m\n) supersymmetry. In this language, the Haldane-Shastry spin chain is of type 
SU(0\2) and the t-J model is of type SU(1\2). While the following results are easily gen- 
eralized to the SU{m\n) case, we are solely concerned with the S77(2|2) extended Hubbard 
model in this letter. 

In addition to the grading of the permutation operator, the coupling parameter Jjj 
determines the nature of this model. The extended Hubbard model studied by Efiler, Ko- 
repin and Schoutens || corresponds to nearest neighbor interactions, J^- = 8iu-j\. We are 
interested instead in the long range version of this model where Jij = l/d 2 (i — j), and 
d(n) = - sin(7m/L). 

Since the Hamiltonian given by Eq. ([!]) acts only as permutations, the total numbers 
of each species are all conserved quantities. We denote the number of localons, unpaired 
up-spin electrons, unpaired down-spin electrons and holes by 

L 

i=l 
L 

i=l 
L 

i=l 
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Q I. - M ■ M - A. (3) 

The total number of unpaired electrons may then be written as M = Mj + Mj_. In the 
following, we work in the Hilbert space of fixed A, M^, Mj_, and Q. The system reduces to 
the S77(l|2) long range t-J model when A = and the S77(0|2) Haldane-Shastry spin chain 
in the limiting case when Q = in addition to A = 0. 

In order to construct the ground state and excitations, we write the states of the system 
in terms of excitations from the fully polarized up-spin state, \F) = niLi c ]f|0)- Thus a 
general state vector is given by 

A Mi Q 

|0) = £ <t>(x,y,z)Ucl d l[ cl^Rc JF). (4) 

x,y,z i=l a =l 1=1 

The amplitude <fi is antisymmetric in the positions of the localons, x = (xi, X2, ■ • ■ , x^), sym- 
metric in the positions of the unpaired down-spins, y — (y\, yi, ■ • ■ , J/m,), and antisymmetric 
in the positions of the the holes, z = (zi, Z2, ■ ■ ■ , zq). When the state vector is written in 
the above manner, the singly-occupied up-spin states correspond to the sites unoccupied by 
the localons, holes, and singly-occupied down-spin electrons. The fact that all four species 
completely fill the lattice has proven useful in allowing alternate equivalent representations 



in terms of excitations from uniform backgrounds of any of the three other species [14 



Assuming that S z = (Mj — M|)/2 > 0, we propose the following Jastrow wavefunction 
= e ^(E^+E^)+^E a ^] x J] d 2 (y a - y p ) x 

l<a</3<M l 

Mi A Q 

x 11(11 d (y<* - x i) • n d (y<* ~ z i)) x II d(xi-Xj)x 

a=l i=l 1=1 l<i<j<A 

X Yl d ( Z l ~ Z n) X II d ( X i ~ Z l)i ( 5 ) 
l<2<n<Q l,i 

where the quantum numbers Jh and J s are either integers or half-integers as appropriate 
so that the wavefunction is single-valued under any of the individual translation operations 
Xi — > Xi + L, y a — > y a + L or z\ — > Zi+L. The Jastrow wavefunction is constructed to satisfy 
the appropriate symmetries upon interchange of any pair of sites. 
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Inspection of the wavefunction indicates that the localons and the holes play the same 
role. This may be understood from the SU(2) symmetry relating the two bosonic states. 
Furthermore, Sutherland has shown that the ground state energy for a general SU(m\n) 



permutation model coincides with that of the equivalent SU(l\n) model P0 |. We therefore 
can take identical localon and hole momenta, Jh, and can then compute the effects of the 
Hamiltonian acting on the wavefunction as if one were dealing with the S77(l|2) supersym- 
metric t-J model of 1/r 2 hopping and exchange. The number of "holes" Q in the effective 
t-J model is given by the total number of bosons, Q = Q + A. 

Using this analogy to the t-J model, we have computed the action of the Hamiltonian 
H on the Jastrow wavefunction defined by Eq. (J5|). We find that the wavefunction is in 
fact an exact eigenstate of the system provided the quantum numbers J s and Jh satisfy the 
following constraints: 

\J h -L/2\ < (M T + l)/2 
\J s -L/2\ <l + S z /2 

\J S -Jh\ <(M l + l)/2. (6) 

Under these conditions, the energy of this state, as a function of the quantum numbers Jh 
and J,, is given by 

E/(n 2 /L 2 ) = h(L 2 - 1) + l -Q{M l + Q)(2M l - Q) 
+M l [^(Mf-l)-2J s (L-J s )] 
+ Q [\{Q 2 -l)-2J h (L-J h )] 

+ 2Q(J S -J h ) 2 . (7) 

We recall that the above eigenenergies have been calculated in a sector of fixed M^, Mj_, 
Q, and A. The ground state of the Hamiltonian in the Hilbert space of fixed Q and A is 
then given by the Jastrow wavefunction in the S z = or 1/2 sector with J s and Jh both as 
close to L/2 as possible. Since the ground state energy is identical among all the Q + A + 1 
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states with total boson number Q, the ground state falls in either the (1,Q + A + 1) or 
the (2,Q + A+ 1) representation of SU(2) x £17(2) C S77(2|2) for S z = or S z = 1/2 
respectively. As a consequence, the lowest energy state in the complete Hilbert space is 
given by Q = L and is L + 1 fold degenerate. The corresponding state vectors are given 

by \G) = E{ x y nil 4 ( t4,i 

|0), where A = 0, 1, 2, . . . , L. Any of these ground states with 
A > is superconducting and exhibits off-diagonal long range order. The state vectors are 
very similar to the 77-pairing eigenstates for the conventional Hubbard model constructed by 
Yang H], except that the phase of each localon is 1. 

The Jastrow wavefunctions corresponding to other values of J s and Jh will be excited 
states of the system. One may compute the corresponding spin and charge excitation energies 
relative to the ground state. In particular, for the Jastrow wavefunctions written in this 
way, various correlation functions, such as the correlation functions between the localons, 
the holes, and the singly-occupied down spins, can be computed in a compact form following 
the recent work of Forrester [^TJ . 

For the extended Hubbard model, the above wavefunctions do not span the complete 
energy spectrum of the system. This comes as no surprise since the Jastrow product only 
describes uniform excitations of the various particles. Nevertheless, since we have shown the 
ground state of this model to be a product of two-body functions, we are led to a solution in 
the form of an asymptotic Bethe ansatz (ABA) [(Tjjlll . With two bosonic and two fermionic 
species to consider, the nesting of the Bethe ansatz may be performed in various manners 
. We choose the BBFF grading and start from the fully polarized background \F). 
Introducing three sets of pseudomomenta, {p a }, {if} an d {?,• }, the ABA equations are 
of the form 

My+Q+A Q+A 

p a L = 2ixJ a + 7T sgn(p a - p p ) - it sg n (P« - <?« (1) ) 

13=1 1=1 
Mi+Q+A a 

2tt// 1) = 7T s g n (<?i (1) - Pa) ~ 7T ^ sgn(^ (1) - qf ') 

a=l i=l 
Q+A A 

27T/j 2) = tt J2 sgn(gf } - - vr £ sgn(g| 2) - qf). (8) 
1=1 3=1 
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The resulting energy and momentum depend only on the {p a } and are given by 



1 Mi+Q+A 
Z a=l 



P= E (P«-*0, (9) 

a=l 

where the reference energy of the system is f £(1 — jjs), and the reference momentum of the 
system is (L — l)n. 

The (individually non-overlapping) quantum numbers J a , ij 1 and 1$ are either integral 
or half-integral and characterize the ABA state. They are restricted to lie in the range 

\J a \ < {L-M l -l)/2 

l4 2) |<(Q-l)/2- (10) 

As a result, in this BBFF picture, the ABA equations only hold in the sectors of the Hilbert 
space satisfying > > A + 1 and Q > A. Empirical observations on small systems 
(L < 10) indicate that the ABA generates highest-weight states of the S'?7(2|2) superalgebra. 
This fact may be used to construct entire states of the extended Hubbard model. However, 
simple counting rules based on spinon degrees of freedom (zeros in J a ) 0,0] no longer 



generate the entire states when additional species are present. In order to account for the 
complete states, it is important to modify the ABA equations by taking the "squeezed" - 
string picture into account ||23[| . Since the resulting Ha-Haldane equations are similar to 
the exact Bethe ansatz solution of the nearest neighbor model 0, the counting of states 
proceeds in a similar manner p^j25[| . 

In both the Haldane-Shastry model and the supersymmetric t-J model, the ABA is 
known to give exact results, even in the non-asymptotic limit |7| JTl|ri^| . While it would be a 
more complicated procedure to prove the same for Eqs. (|8|), the integrability of this model 
as well as its similarities to the t-J model should enable such a proof. 

In summary, we have provided the ground state and uniform excitations of the long range 
extended Hubbard model. These wavefunctions are in the form of a Jastrow product, for 



which various correlators are already available. This system of strongly correlated electrons 
is completely integrable with an infinite set of mutually commuting constants of motion. 
This fact enables us to determine the full energy spectrum in terms of an asymptotic Bethe 
ansatz and its Ha-Haldane modification. As a result, it should be possible to analytically 
study the thermodynamics and elementary excitations of this model. Since we only needed 
to know the general properties of the graded permutation operator ILj, our results are easily 
generalized to the arbitrary long range SU(m\n) model. 

This work was supported in part by the Swiss National Science Foundation and by the 
U. S. Department of Energy under grant no. DOE-91ER40651-TASKB. 
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